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Usually alternative theories of gravity imply deviations from the well-known Kerr space-time, a 
model of an isolated black hole in General Relativity. In the dominant order, the deformed Kerr 
metric, free of closed time-like curves outside the event horizon, has been suggested recently by Jo- 
hannsen and Psaltis. It has a single deformation parameter which is not constrained by the current 
observations, allowing, thereby, for a kind of unified and simple phenomenological description of 
black holes in various theories of gravity. Here we consider a number of classical and quantum phe- 
nomena of radiation in the vicinity of such deformed Schwarzschild-like and Kerr-like black holes: 
spiralling of particles into black holes, decay of fields propagating in the black hole's background. 
Hawking radiation. In particular, we calculate some quantitative characteristics of the above phe- 
nomena, such as the binding energy of particles, quasinormal modes, late-time tails of fields of 
various spin, intensity of Hawking radiation. The binding energy released when a particle goes over 
from a given stable orbit in the equatorial plane to the innermost stable one is calculated for such 
non-Kerr black holes. Due to inseparability of the wave equations in the general case, the pertur- 
bations and stability of scalar, Dirac, and electromagnetic fields are analyzed for vanishing rotation 
only. The dependence of the radiation phenomena on the deformation parameter is discussed. 

PACS numbers: 04.30. Nk,04.70.Bw 



I. INTRODUCTION 

Rotating black holes in 3-f 1 General Relativity are de- 
scribed by the well-known Kerr solution, which is singled 
out by the uniqueness theorem: Kerr space-time is the 
only stationary, axisymmetric, asymptotically flat solu- 
tion to the Einstein equations in vacuum, which possesses 
an event horizon and free of closed time-like curves out- 
side of it. Thus, only mass M and angular momentum 
L completely determine the properties of black holes in 
General Relativity. Therefore, Kerr space-time has been 
extensively studied since its discovery in the context of 
various potentially observable phenomena, such as radia- 
tion of gravitational waves, accretion of matter, gravita- 
tional lensing, and also in the scope of theoretical ques- 
tions related to evolution of black holes, gravitational 
collapse, Hawking radiation and others [ij. 

At the same time, Einstein's General Relativity, being 
apparently the simplest geometric theory of gravitation 
among others, raises a number of fundamental questions 
concerning the nature of singularities, dark energy and 
dark matter, quantization of gravitational interactions. 
Attempts to resolve these questions lead theorists to a 
number of alternative theories of gravity, which give the 
same observational consequences at the modern level of 
experiments, because Einstein's gravity, as well as a num- 
ber of its alternatives, has been tested only in the weak 
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field regime, i.e. in the post-Newtonian approximation. 
A future experimental test of the strong field regime may 
be given by gravitational antennas, which could detect 
gravitational waves from neutron stars and black holes 
through observation of the quasinormal modes of these 
compact objects. Recently a number of works have been 
devoted to the study of potentially observable processes 
around various non-Kerr black holes pl-^. Therefore, a 
unified description of analogues of the Schwarzschild and 
Kerr solutions in various theories of gravity would be 
most useful for testing the alternatives. Such a model 
was suggested by Johannsen and Psaltis [Tj who con- 
sidered deviations from the Schwarzschild and Kerr so- 
lutions and found a regular outside the event horizon 
space-time, which reduces to the Kerr one, when the de- 
formation parameters vanish. Current observations re- 
strain values of the deformation parameters, so that in 
the dominant order only one deformation parameter e re- 
mains which leads either to more prolate (e > 0) or more 
oblate (e < 0) shape of the horizon than the one of Kerr. 
Thus, the Johannsen-Psaltis space-time is described by 
three parameters: mass M, angular momentum per unit 
of mass a and deformation e. 

Recently the Johannsen-Psaltis background has been 
studied in a number of papers [sl-fl^ with a special em- 
phasis on gravitational lensing and Penrose processes of 
energy extraction around such black holes. Here, our aim 
is to study various processes of radiation (both classical 
and quantum) around such a generic black hole's model. 
A classical radiation of a field near a black hole consists 
of three stages, at each of which one of the following pro- 
cesses is dominating: 1) initial outburst, which crucially 
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depends on the initial perturbation, 2) the intermedi- 
ate damped oscillations represented by complex frequen- 
cies, termed quasinormal modes, and 3) the asymptotic 
tails at late times. Therefore, here we study both pro- 
cesses which are independent on the initial conditions of 
the perturbation, the quasinormal modes and late time 
tales, though, in the limit of vanishing rotation, due to 
inseparability of variables in field equations for the most 
general case. The Hawking radiation, for estimation of 
which one requires the classical reflection index of fields, 
is also considered for the non-Schwarzschild case. 

Another important characteristic which we consider is 
the binding energy of a particle moving from a given 
equatorial orbit to the innermost stable one. The binding 
energy allows one to learn how much energy the matter 
(for example an accretion disk) will release before plung- 
ing into the black hole. Here we compute the binding en- 
ergy for arbitrary rotation and deformation parameters a 
and s. We have found that the binding energy is increas- 
ing if the shape of the horizon is more prolate than the 
one of Kerr solution, that is, when e < 0. The Hawking 
radiation, on the contrary, is enhanced for larger values 
of e. Real oscillation frequencies of quasinormal ringing 
are increasing when the deformation parameter e grows. 

The paper is organized as follows: In Sec H we give 



some basic properties of the Johannsen-Psaltis black hole 
which will be explored in Sec HI, devoted to deducing the 
wave equations and stability of scalar, electromagnetic, 
and Dirac fields in the spherically symmetric background. 
Sec. IV describes numerical methods used in the paper. 
In order to find quasinormal modes we used Frobenius 
expansion fis]. 16j , WKB method , and time-domain 
integration [18 [, the last one was also used for getting 
late-time tails. Sec V summarize the obtained results on 
quasinormal modes and late-time tails, while in Sec VI 
the intensity of Hawking radiation is computed for fields 
of various spin. 

Note that throughout the paper we use geometrical 
units G = c = 1 and fix the scale symmetry in all 
the numerical calculations such that the black hole mass 
M = 1/2. 



II. THE NON-SCHWARZSCHILD AND 
NON-KERR BACKGROUNDS 



In Boyer-Lindquist coordinates the Johannsen-Psaltis 
metric [tI can be written as 



[l + h{r,e)] [ 1 



2Mr 



AaMr sin^ 9 
+ E 



where — + a?' cos^ 6*, M is the black hole mass, a is 
the rotation parameter, and the deformation parameter 
e comes from the following general expansion 



M^,^)-f:(^2.+e..+,§:)(^)'. (2) 
fc=o ^ / V / 

The current Lunar Laser Ranging experiment con- 
strains the post-Newtonian parameters of the metric, so 



[1 + h{r, e)]dtd(j> - 
+ h{r, 9) 



S[l + fe(r, 9)] 
A + a2 sin^ 9h{r, 9) 
a'^{ll + 2Mr) 9^ 



dr^ 



r 



that, in the dominant order, one can take 



h{r,9) = e- 



(1) 



(3) 



Various astrophysical constraints on e can be also found 
in 0,113. 

The inverse metric has the following form: 



/ r^{r'^W{2M+r)[eM'^+r^)) 
I (eA/a+ra)((r-2M)r4+a^(eAf-*+r-i)) ^ 






I 2aMr^ 

\ {r-2M)r^+a'-'{sM'-i+r'-') 



2aMr^ \ 

{r-2M)r*+a'^{sM'-*+r-') ' 



(4) 



r^(r-2M) i 



The event horizon is located at the corresponding root of rameter e, once positive (negative) leads to a more pro- 
the equation g^^ — §1194,4, ~ 0, and the deformation pa- 
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late (oblate) object than the Kerr black hole. As ggg > 
and grr > 0, g^^ > 0, this spacetime is free of closed 
timelike curves. For arbitrary e < the space-time has 
a regular event horizon and represents therefore a black 
hole. On the contrary, e > constrains the maximal 
value of the rotation parameter a (see Fig. 2 in 7] ) , which 
becomes smaller than the one of Kerr solution. 



III. BINDING ENERGY 

In order to learn how much energy could be extracted 
from a particle slowly spiralling toward the black hole's 
horizon, it is necessary to learn what the energy of the 
particle at the innermost stable circular orbit is. For this 
purpose it is useful to introduce a function termed the 
binding energy, which is the amount of energy released by 
the particle going over from a given stable circular orbit 
located at Tc to the innermost stable orbit r^^. Thus, the 
binding energy in percent is 



Binding Energy = 100 x 



E{rc) 



(5) 



Frequently, instead of some given circular orbit Vc a par- 
ticle coming from infinity is considered, so that Tc = oo is 
taken in the above definition of the binding energy. The 
binding energy for the Schwarzschild black hole is about 
5.7 percent and for the extremal Kerr black hole can be 
as large as 42 percent. This is much larger than, for in- 
stance, the binding energy at the nuclear fusion, which is 
about 0.7 percent of the rest mass energy. Recently the 
binding energy for charged particles around a black hole 
deformed by the tidal force and magnetic field has been 
found in iiHI. 

For our case of a stationary, axisymmetric metric, one 
has three integrals of the motion. 



P 



which are particle's rest mass /i, energy E — pt, 
angular momentum L = — p^. Then, the equation 



PaP 



(6) 



and 



(7) 



leads in the equatorial plane to the following relation: 



1 

Qrr 



[g'^E^ - 2g"^EL + g**L^ - fj.^) . (8) 



The energy E and momentum L of a particle at a circular 
orbit r can be determined numerically from the require- 
ments: 



dr ^ rf^r ^ 
dr ' dr^ 



(9) 



Then, the numerically obtained E{r,a,e)/fi, L{r,a,e)/fi 
allow one to use formula ([S]). 



TABLE I: Binding energy Be for various values of the de- 
formation parameter e and rotation a. Left columns are for 
particles co-rotating with a black hole, right columns are for 
counter rotating particles. The star means that the extremal 
values a — 0.5 for co-rotating particles was reached in com- 
putation only asymptotically, as the innermost stable orbit 
approaches the horizon in this case. 





CO- rotating 


counter rotating 


e 


a 




E{ns) 


Be 






Be 


-5 


0.1 


3.26251 


0.947203 


5.27972 


3.78356 


0.954798 


4.52022 


-5 


0.2 


3.01354 


0.942488 


5.75122 


4.04775 


0.957840 


4.21596 


-5 


0.3 


2.78166 


0.937117 


6.28830 


4.31171 


0.960488 


3.95116 


-5 


0.4 


2.57733 


0.931213 


6.87875 


4.57452 


0.962809 


3.71907 


-5 


0.5* 


2.41140 


0.925130 


7.48702 


4.83569 


0.964858 


3.51420 


-1 


0.1 


2.79356 


0.938409 


6.15909 


3.41641 


0.949953 


5.00469 


-1 


0.2 


2.46574 


0.929858 


7.01424 


3.71439 


0.954065 


4.59345 


-1 


0.3 


2.12623 


0.917971 


8.20290 


4.00518 


0.957474 


4.25260 


-1 


0.4 


1.78144 


0.900507 


9.94928 


4.28970 


0.960353 


3.96472 


-1 


0.5* 


1.46393 


0.874561 


12.5439 


4.56904 


0.962822 


3.71779 


1 


0.1 


2.53274 


0.931827 


6.81727 


3.22108 


0.946834 


5.31664 


1 


0.2 


2.14187 


0.918816 


8.11845 


3.53857 


0.951732 


4.82680 


1 


0.3 


1.67902 


0.895150 


10.4850 


3.84436 


0.955666 


4.43339 


5 


0.1 


2.04540 


0.911207 


8.87926 


2.82637 


0.938547 


6.14526 


5 


0.2 


1.59086 


0.877896 


12.2104 


3.17842 


0.945817 


5.41832 


5 


0.25 


1.30398 


0.841900 


15.8100 


3.34747 


0.948706 


5.12938 



On Fig. [T] one can see that the more prolate (or the 
less oblate) shape of the horizon is, the larger the bind- 
ing energy is. This is true for any values of the rotation 
parameter a, as it is shown on table ID An increase in the 
binding energy is slightly larger for highly rotating black 
holes than for slowly rotating ones and can be about 100 
percent for relatively small deformations. Indeed, e = 5 
means deformation of order eM^/r^ ^ 5/8 and corre- 
sponds to approximately 2.5 times increase (see table H]) 
in the binding energy. 

We know that the binding energy for the extremal Kerr 
black hole equals 3.8 percent for counter-rotating parti- 
cles and 42 percent for co-rotating ones. In concordance 
with this we obtain in the limit of small e for the non-Kerr 
black hole: for a = 0.4999999, e = -lO'^, Eb = 41.46 
percent for co-rotating particle and Ei, = 3.78 percent for 
counter rotating one. Thus, in the limit of vanishing de- 
formation parameter the binding energy approaches its 
Kerr values. 



IV. DEDUCTION OF THE MASTER WAVE 

EQUATIONS FOR SCALAR, 
ELECTROMAGNETIC, AND DIRAC FIELDS 

In the most general case we were unable to decouple 
variables in the field equation, and had to be limited 
by the non-rotating case. Although scalar and electro- 
magnetic wave equations can be easily derived for this 



4 



binding Energy 
12 h 



10 
8 
6 
4 
2 



Binding Energy 




20 40 



60 



100 




60 80 100 



FIG. 1: 1) Left panel: Binding energy as a function of Tc for e = 5 (top), 0.01, and —5 (bottom), a = 0.2 for particles co-rotating 
with a black hole at a circular orbit rc- 2) Right panel: Binding energy for e = 5 (top), 0.01, and —5 (bottom), a = 0.2 for 
counter-rotating particles at a circular orbit rc. 



case, the Dirac field requires sorae algebraic calculations. 
Therefore, we deduce here the wave equations for spin 0, 
±1/2 and ±1 fields in the general spherically symmetric 
background using the tetrad formalism. 



A. Spherically symmetric backgrounds and the 
tetrad formalism 



A spherically symmetric background can be described 
by the following line element 

ds"^ = g^.^dx^'dx" = Adf - —dr^ - r^{de^ + sin^ 9dip'^), 

(10) 

where the metric coefficients A and B depend in general 
case on both t and r coordinates. 

The scalar field satisfies the Klein-Gordon equation. 



The massless Dirac and Maxwell equations are 



(11) 



(12) 



\/2Vbb'P^ = 0, 
\/2Vbs'Q^ = 0, 
^^■M = {F^,, = d^A, - a,A^), (13) 



where and Vbs' denote covariant differentiation. 

Following [IJ we rewrite (fT2)) and ([T3| in a tetrad sys- 
tem of null- vectors: 



h = ^{A,B,0,0}, m, 
B , _ 



= ^{0,0, l,isin6'}, 



V2 
r 

71 



{0,0, l,-isin6'}, 



(14) 



that satisfy the following relations 

l^jl^ — n^ji^ — m^jn^ = rrifjn'^ = 0, 
l^mf^ — Z^m^ — n^m^ — rif^m^ — 0, 

Z^rt^ = 1, mfjn'' = -1, 
l^n^ + rifj^ly - m^Jfli, - m^m,, = g^j_^. 

The spin coefficients are 

K = If^.^m'^r =0, TT = -n^.^rn'^r = 0, 



B-A' 
4B ' 



A 



A = —nfj_-^m^m'^ = 



2rB' 

a = -{l^i-un^rrf - m^-iM^^m!') = — (15) 



2\/2r' 



a = l^j^-jjin^m^ = 0, 
1 



/i = —rifj^.um^m 



1 

Vb' 



cot 6' 
2V2r' 



-n^;,j7n^n" = 0, r = l^j^.^i^m^n^ = 



2^2 2AB ' 

where dot and prime denote derivatives with respect to 
t and r. 

The intrinsic derivatives read 



A = 



(5 = ;^m^ 



(9 



(16) 



<5 = ;^m'' = - 



x/2r96l v^rsin^cl^' 
19 id 



V2rde ^/2rsmed(p' 
In terms of these variables the massless Dirac equations 
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can be reduced to [24 



Substituting the following ansatz 



{D + e-p)P° + (6 + 7r-a)P'^ = 
(A + /i - 7)pi + ((5 + /3 - T)P0 = 
(D + e - p)q" + (5 + TT - S)Q^ = 
(A+7I-7)Q' + (5 + ^-r)Q" = 

The Maxwell equations take the form [23j : 



D^i-S^o = (vr - 2a)$o + 2p<i>i - K<i>2, 

£)$2-?$i = -A$o + 27r$i + (p- 2e)>I>2, 

5$i - A$o = (^t - 27)<I>o + 2r$i - a-$2, 

(5$2 - A$i = -j/<I>o + 2^$i + (r - 2/3)$2, 



(17) 



(18) 



where $o = F^J^'m" , $i = \F^,^{l^'n'' + 777^771"), and 



+ l/2>£™(^^,¥')^-'/'^(i,0"'/', 
-l/2^£m( 



$2 - F^^mi'n\ 



R+i/2{t,r) 

R-iMt, r) .y^Yiraie, v)r-^/^A{t, rf/^ 

R-i/2it, r) +i/2r,„(0, ^)r-^/^A{t, ry/\ 
R-i{t,r) ^,Ye^{e,ip)r-^A{t,ry^^, 
R+i{t, r) +iYem{0, v)r~^A{t, r^^/^ 

(19) 

into pill7ll8l) (here, sYem{9,(f>) are the spin- weighted 
spherical harmonics with £ — \s\, \s\ + 1, |s| + 2, . . . and 
m = —£, ...£—!,£), we obtain the following second 

order equation for Rs(t,r): 



pO 

$0 
$2 



A d 



(AdRs 


h 


AB - 2AB dRs 


\B dr 




' A^ 


2AA' 


2A^B' AA'B- 


^252 + 


rB^ 


rB^ 2B3 


A^ 


AA' 


A'^ 2A^B' 




rB^ 


^ 452 rB^ 




(20) 



Rs = 0, 



where A = (^ + |s|)(^ - |s| + 1). 

Although we have derived equation (PU]) for s — 
0,±l/2, ±1, we will use pO| only for scalar and Dirac 
fields, while for electromagnetic perturbations we shall 
use the following simpler equation, 



A d A dR±i \ A d B dR±i \ AXR±i 



B dr \ B dr 



Bdt\ A dt 



= 0, 



(21) 

which can be derived from either the first and third or the 
second and forth equations of (jlSp by the substitution 



$1 = R±i{t,r) ±iYi^{e,ip)r- 



(22) 



B. Linear perturbation equations and stability of 
fields in the non-Schwarzschild black hole 
background 

The Schwarzschild-like metric takes the form 

ds^ = (l + hir)) (^fir)dt^ - -r2(d02+sin2 0V), 

(23) 

2M 

where /(r) = 1 , h{r) — e — ^. 



Making use of the relations A ^ {1 + h{r))f{r), B = 
1 + h{r), and Rs{t,r) = e-'"*^'s(r) in wc find the 

wave- like equation 



drl 



(24) 



where 



dr 



is the tortoise coordinate. The effective potential Vg for 
scalar, Dirac, and Maxwell fields has the following form 



Vo = k\ + 



f{r)f{r) 



V+i - fcA, 



(25) 



dk 



— iuj dk 1 d / 1 dk \ 1 
^+5= ^2^d^^Ad^[kd^)^16k^\dr^ 



1 + h(r) 

where k = f{r). Note that for h{r) =0 the po- 

tential V_^_l coincides with (2.17) of [25|. We will consider 
only s > because of the symmetry between opposite 
polarizations in a spherically symmetric background. 

The effective potential for scalar and Maxwell fields 
has the form of the potential barrier that vanishes at the 



6 



Voir) 




FIG. 2: The effective potentials {£ — 1) for the scalar (left panel) and Maxwell (right panel) fields for e — —4 (red, lower), 
£ = (green) and £ = 4 (blue, upper). 



event horizon and infinity and is positive definite every- 
where outside the event horizon (see fig. [2]). Therefore, 
the solutions to the wave-hke equations are eigenvectors 
of a positive self-adjoint operator in the Hilbert space of 
square integrable functions, implying that there are no 
growing modes in the quasinormal spectrum. 

The effective potential for Dirac fields is complex due 
to its dependence on w. However, similarly to the 
Schwarzschild black hole, one can introduce an equiva- 
lent wave-like equation with an effective potential that 
contains a square root. Following [2^, we define a new 
function Z, such that 



ILU ] Z + . 



(26) 



Substituting ^^+1 into the wave-like equation ([24| we find 
that the new function Z satisfies the wave-like equation 
with a real effective potential 



fz_ 

drl 



- Xk 



dVXk 
dri. 



Z = 0. 



(27) 



The corresponding differential operator is also positive 
self-adjoint one what can be shown with the help of S- 
deformation [23| by taking S — Vxk. Let us note, that 
we shall not use Eq. (I?fl) for finding quasinormal modes, 
because the Frobenius method cannot be easily applied 
to potentials with square roots [1^ [11] . 

Thus, the stability of massless fields in the nonrotating 
Johannsen-Psaltis background is straightforward ^29] and 
basically identical to the Schwarzschild case. 

The potential barrier (fig. [3]) is growing when the de- 
formation £ is increasing. As we shall see later this fea- 
ture of the effective potential will stipulate the character 
of dependence of the quasinormal modes and scattering 
amplitudes on e. 



V. NUMERICAL METHODS 

In this section we briefly discuss the numerical meth- 
ods used for calculations of the quasinormal modes, the 
transmission/reflection coefficients, and late-time tails. 



A. Quasinormal modes and the Frobenius method 

In order to calculate quasinormal modes we impose 
the quasinormal mode boundary conditions for the wave 
equation (j24p . that is, we require that there are only 
purely ingoing waves at the black hole's horizon, 

^'(r* — —oo) (X exp(— iwr*), 

and only purely outgoing waves at spatial infinity, i.e. 

^'(r* — !• oo) cx exp(+«Ci;r*). 

Thus, no waves are coming from the horizon or infinity, 
which implies that uj are proper oscillation modes in the 
black hole response to an "instantaneous" perturbation. 
In other words, when the perturbation decays, the source 
of the initial perturbation is not acting anymore. 

Equation ([M)) has an irregular singularity at spatial 
infinity and a regular singular point at the horizon r = 
2M. The appropriate Frobenius series has the form 



2M 



R 



-ib 



Rf 



n=0 



2A/\" 



where R is the closest to the horizon singular point. 
We find that h — 2Muj, c — 2Muj for the scalar and 
Maxwell fields, while for the Dirac field b = 2Muj — i/A 
and c — 2Muj + i/2, because of the w-dependent effective 
potential. 

For the scalar and Maxwell fields the only other singu- 
larity is r = and, therefore, R = 0. For the Dirac field 
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there are additional singular points which are solutions 
to equation 1 + h{r) — 0, so that R — for £ > 0, and 
R = M-^ for £ < 0. 

The coefficients a„ satisfy a recurrence relation which 
can always be reduced to the three-term one through the 
Gaussian eliminations and, finally, one can find the equa- 
tion with an infinite continued fraction (see e.g. |29j] for 
details). This equation has an infinite number of roots u;„ 
corresponding to the quasi- normal frequencies p^ . The 
infinite continued fraction converges very slowly when the 
imaginary part of w is large. In this case, in order to im- 
prove the convergence, one can use the Nollert procedure 

m 

B. Reflection coefficients 

In order to calculate the Hawking emission rates of par- 
ticles, one needs to solve the problem of classical scatter- 
ing and find the gray-body factors of the corresponding 
fields. Such a problem implies classical scattering bound- 
ary conditions for the wave equations obtained in Sec 
III. Thus, at the event horizon the boundary condition 
corresponds to a purely ingoing wave, while at spatial in- 
finity (r — > cx)) one has a sum of the ingoing and outgoing 
waves, 

*(r) ~ Zin exp(-iiLjr*) + Zout exp(za;r^), 

where Zin and Zo^t sue integration constants. Thus, we 
would like to know which portion of particles will be able 
to pass through the barrier of the effective potential. 
Introducing the new function 

and choosing the integration constant as P{2M) — 1, 
we expand Eq. (j24[) near the event horizon and find 
P'(2M), what completely fixes the initial conditions for 
the numerical integration. Then, we integrate Eq. (j24p 
numerically from the event horizon 2M to some distant 
point rf ^ M and find a fit for the numerical solution 
far from the black hole in the following form: 

P{r) = Z.nPUr) + ZoutPoutir), (28) 

where the asymptotic expansions for the corresponding 
functions are found by expanding (1241) at large r as 

P,„(r) = (1 + 4:^-1 +/^l?V-2 + ...), 

Poutir) = e-V- (l + P^r-' + p!,llr-' + ...). 

The fitting procedure allows us to find Zin and Zout- In 
order to check the accuracy of the found coefficients, one 
should increase the internal precision of the numerical 
integration, the value of r^, and the number of terms in 
the series expansion for Pin{r) and Poutir), making sure 



that the values of Zin and Zout do not change within de- 
sired precision. Thus, the above approach can be used for 
finding reflection coefficients of a wide class of compact 
objects with great accuracy (soj . 

C. WKB method 

In order to obtain an analytical expression for the 
quasinormal frequencies in the regime of large multipole 
numbers £ and also for an additional check of the data ob- 
tained by the convergent Frobenius method, we shall use 
the WKB formula of the 6th order beyond the eikonal 
approximation fvA [3l| . The formula has the following 
form: 

!(^^-gA..n+i, ..0,1,2..., (29) 

where the correction terms were obtained in [TtI . Isij 
and depend on higher derivatives of V at its maximum 
with respect to the tortoise coordinate r*, (n labels the 
overtones). The WKB approach was developed by Schutz 
and Will [17] and extended up to the 3rd and 6th orders 
in [31i] . It can be effectively used not only for finding low- 
lying quasinormal modes, that is, oscillations with longer 
lifetime (see for instance [32?! and references therein), but 
also for calculations of the transmission/reflection coeffi- 
cients in various problems [1^ [sj] ■ 

D. Time-domain integration 

In addition, we shall study late-time tails through the 
numerical characteristic integration method that uses the 
light-cone variables u = t — r*,v = t-\-r*. In the charac- 
teristic initial value problem, initial data are specified on 
the two null surfaces u = uo and v = vq. The discretiza- 
tion scheme was suggested in [l^ and used in a number 
of subsequent papers (see, for instance [111), showing an 
excellent concordance with accurate Frobenius data at 
the stage of quasinormal ringing. 

VI. QUASINORMAL MODES AND LATE-TIME 
TAILS 

Quasinormal modes can be written as 

w„ = Re{ujn) + ilm{ujn), (30) 

where n is the overtone number. The real part of the fre- 
quency is the real oscillation frequency, while the imag- 
inary part is proportional to the damping rate of the 
mode. In our designations, once /m(a;„) is negative, the 
mode is damped, while positive imaginary part corre- 
sponds to growing (unstable) modes. As at late time 
the dominant in a signal mode corresponds to n = 0, 
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FIG. 3: Real and imaginary parts of the fundamental quasinormal mode as a function of e for the scalar field 1=1 (blue line, 
left panel top, right panel bottom), Dirac field I = 1/2 (green line, left panel bottom), and Maxwell field 1=1 (red line, right 
panel top) fields. 



we shall consider here only this lowest mode, termed the 
fundamental modes. 

On Fig. |3] the fundamental quasinormal modes (the 
modes with the smallest decay rate) computed by the 
Frobenius method are shown for all the fields under con- 
sideration. A striking feature of the fundamental modes 
are almost linear dependence on e of their real oscillation 
frequencies. Possibly, this is connected with a relative 



smallness of the deformation h for e ranging from —5 to 
5. In the regime of strong deformations, that is certainly 
not justified physically, the dependence on e might be 
not linear anymore. The fundamental modes can be ap- 
proximated by the following analytic expressions for the 
scalar (£ = 1), Dirac {i = 1/2), and Maxwell {I = 1) 
fields respectively: 



uj w 0.58587(1 -I- 0.01756e-f0.00000e2)- 0.19532i(l-f 0.00363£-f0.00045£2) (scalar) 

u « 0.36593(1 + 0.02179e- 0.00012£2)_ 0.19397i(l- 0.00054£ + 0.00134e2) (Dirac) (31) 

uj w 0.49653(1 + 0.02251e + 0.00000e2)- 0.18498i(l-f0.00244e + 0.00101e2) (Maxwell) 



Thus, the real oscillation frequency indeed grows almost 
proportionally to e, as can be concluded from smallness 
of the coefficients in front of . The damping rates, given 
by absolute values of imaginary parts of the frequencies, 
are not monotonic function of e, having a minimum at 
some negative e (different for each type of perturbation) 
corresponding to the least damping. Let us note, that 
as the potential barrier is monotonically growing when e 
increases (fig. [2|), some growth of Re{uj) (certainly, not 
necessarily linear) as a function of e was expected. 

The above numerical data obtained by the Frobenius 
method can be also checked by the time-domain inte- 
gration (see Fig. 2]). By fitting the time-domain pro- 
file we are able to reproduce the fundamental modes ob- 
tained in the frequency domain within any desired ac- 
curacy. The asymptotic late-time tails coincide with the 
Schwarzschild one [36] 



1*1 cx 



t 



-2£-3 
-2£-2 



scalar. Maxwell; 
Dirac; 



t oo. 



(32) 



As time-domain profiles include contributions from all 
modes, this also proves the stability of fields in the in the 
Johannsen-Psaltis background. 

The WKB method is known to be accurate in the 
regime of large multipole numbers and small overtones 
n <ti £. Therefore, for £ 3> 1 the eikonal approximation 
allows us to find simple analytic expressions for QN fre- 
quencies. Expanding the location of the maximum of the 
effective potential rg in powers of 1 /£ and e and making 
use of this expansion in WKB formula (1^^ at the first 
WKB order and lowest orders of e, we have found the 
following expression: 



where 



(33) 



At large i the data given by the above ekinal formula are 
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|R„| l^'ml |K,| 




20 50 100 200 500 20 50 100 200 500 20 50 100 200 



FIG. 4: Time-domain profiles from left to right: for the scalar (i = 1), Dirac {i = 1/2), and Maxwell {£ — 1) fields for e = —4 
(red), e = (green) and e = 4 (blue). Late-time tails decay law does not depend on e. 



in excellent concordance with the numerical one. It is 
well known, that the eikonal formula works well already 
for moderate values of £ = 2,3,4, {n = 0). Comparing 
the analytical expansion ([55]) with (PT|) one can observe 
that the eikonal formula gives quite a good approxima- 
tion for the the real part of w as a function of e already 
for ^ = 1 = 1/2 for the Dirac field). 

It is well known that in the eikonal approximation 
quasinormal frequencies of fields minimally coupled to 
gravity, that is, propagating in the black hole back- 
ground, are determined by the centrifugal term of the 
effective potential, which is independent on the spin of 
the field. Therefore, as a rule, the eikonal expression for 
id (j33p does not depend on spin as well. Thus, we could 
speculate that there may be a kind universal dependence 
on the deformation parameter for minimally coupled to 
gravity fields of any spin, when h is small: the imagi- 
nary part of the fundamental quasinormal mode remains 
relatively close to that of the Schwarzschild black hole, 
while the real part grows proportionally to s. This might 
be true also for gravitational perturbations, if the lowest 
dynamical mode ^ = 2, as it takes place in General Rel- 
ativity. However, if other fields are coupled to gravity, 
such as the dilaton, a spherically symmetric dynamical 
mode £ = appears in the spectrum, which evidently 
cannot be described by the approximation of geometrical 
optics. 



VII. HAWKING RADIATION 

Here we assume that the black hole is in the thermal 
equilibrium with its surroundings in the following sense: 
It is supposed that the black hole temperature does not 
change between the emission of two consequent particles, 
or, in other words, the system forms the canonical ensem- 
ble. 

Once the coefficients Zi„ and Zout are calculated with 
the help of the shooting method related in the Sec V. B, 
one can find the absorption probability 



P — 1 — \Zout/Zin\^ 



(34) 



TABLE II: Energy-emission rates due to radiation of massless 
fields. 



£ 


scalar 


Dirac 


Maxwell 


-4 


0.000312 


0.000213 


0.000102 


-2 


0.000304 


0.000187 


0.000083 





0.000298 


0.000163 


0.000067 


2 


0.000292 


0.000143 


0.000055 


4 


0.000288 


0.000125 


0.000044 



absorption probability for massless Dirac particles we fol- 
low 13711: 



(35) 



In the semi-classical approximation the asymptotically 
flat black hole creates and radiates particles with thermal 
spectrum [38l |. The energy-emission rate is proportional 
to the absorption probability lAi]^ (see also [3§|): 



dM f dio ^ 



{21+1)\A,\''lo 



27r cxp{u;/Th) ± 1 



(36) 



for scalar and Maxwell fields. In order to calculate the 



for particles of each spirality, where Th = 1 /8ttM is the 
Hawking temperature. The infinite sum in p6p converges 
quickly, so that in practice we need to sum over few lowest 
multipoles £ (note that for fermions we sum over half- 
integers £ and take sign in the denominator). 

The results of computations for the energy emission 
rates are shown on Fig. [5] The larger the absolute value 
of e is, the more suppressed emission rates are. On ta- 
ble in] one can see that the deviation from Schwarzschild 
values are larger for fields of higher spin. This can be 
easily understood from the expression for the effective 
potentials of scalar, Maxwell, and Dirac fields Since 
A always appears with the factor k, which contains e, the 
larger multipole number £ is, the bigger deviation arises 
in the potential due to e. At the same time the effec- 
tive potential for the scalar field does not depend at all 
on s when ^ = 0, which corresponds to the dominant 
contribution in the emission rates. 
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VIII. CONCLUSIONS 

Investigation of potentially observable phenomena 
around black holes in a large number of alternative the- 
ories of gravity would be a difhcult, if not never end- 
ing, task, since even solutions for black hole metric are 
not known in many of such theories. A kind of phe- 
nomenologically unified description of black holes by the 
Johannsen-Psaltis space-time gives us the generic model 
which can be studied with simple numerical tools. Here 
we have analyzed a number of radiation processes in the 
vicinity of the Johannsen-Psaltis space-time through con- 
sideration of their essential characteristics: 1) the binding 
energy of a particle spiralling into the black on an equato- 
rial orbit, 2) proper (quasinormal) modes, late-time tails 
and stability of scalar, Dirac, and Maxwell fields, 3) in- 
tensity of the Hawking radiation of these fields. 

We have shown that the binding energy is increasing 
for more prolate than Kerr configurations, what corre- 
sponds to positive deformation parameter e. In a similar 
manner, real oscillation frequencies of quasinormal modes 
are larger for positive e, though, at asymptotically late 
time the power-law tails coincide with the Schwarzschild 
ones. We were able to find analytic expressions for the 
frequency in two regimes: 1) for the fundamental modes 



and small and moderate e and 2) in the eikonal regime 
£ ^ 1. We have found that the real oscillation frequen- 
cies are approximately linear in e, for small and moderate 
values of e. Stability of the fields under consideration has 
been proved. The intensity of the Hawking radiation is 
shown to be suppressed for positive e and enhanced for 
negative ones. 

Our work could be complemented by the analysis of 
classical and quantum radiation processes for massive 
fields. A drawback of our approach is in the impossibility 
to decouple variables in the field equations for a non-zero 
rotation parameter, so that quasinormal modes and the 
Hawking radiation were analyzed only for non-rotating 
black holes. At the same time, strong dependence of 
the QNMs and intensity of Hawking radiation on the de- 
formation parameter e leaves us hope that the obtained 
results will remain qualitatively the same at least for slow 
rotations. 
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